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CHAPTER / 


The Hoyle-Narlikar Theory 


of Gravitation 


E LncroauctiLon 


The success of Maxwell's equations has led to 
electrodynamics being normally formulated in terms of fields 
that have degrees of freedom independent of the particles in 
them. However, Gauss suggested that an action-at-a-distance 
taeory in which the action travelled at a finite velocity 
might be possible. This idea was developed by Wheeler and 
E se i 5 ; TENC 
Feynman (1,2) who derived their theory from an action-principle 
that involved only direct interactions between pairs of part- 
icles. «a feature of this theory was that the 'pseudo'-fields 
introduced are the half-retarded plus half-advenced fields 
claculated from the world-lines of the particles. However, 

1 | | | m B) 
Wheeler and Feynman, and, in a different way, Hogarth 
were able to show that, provided certain cosmological 
conditions were satisfied, these fields could combine to 

i 1 " 4 5 
cive the observed field. Hoyle and Narlikar (4) extended tne 
tueory to general space-times and obtained similar theories 
: ip ذم ابم(‎ (5) | TUN : (6) 
for their 'C'-field and for the gravitational field è 


It is with these theories that this chapter is concerned. 


It will be shown that in an exp.ndinz universe the 


advanced fields are infinite,and the retarded fields finite. 
Tais is because, unlike electric charges, all masses have the 


same Sign. 


2. he Boundary Condition 
Hoyle and Narlikar derive their theory from the 


action: 


A = 2 2. ( (a,b) dadb, 
يه‎ fb 
where the integration is over the world-lines of particles 
Q, NTC In this expression ^ is a Green function 


that satisfies the wave equation: 


C(x XK"). qlr GROKX) < d 
-9 


where q is the determinant of 9: . oince the double sum 
in tne action A is symmetrical between all pairs of 
particles @,b , only that part of G (a,b) that is 
symmetrical between « and b will contribute to the action 
i.e. tne action can be written 
Ls r ١ ‘ ; 
/ = £2 (IC (a,b) a, d b 
where Gla b) $C (a,b) + 40 ), o) 
Ue GO. - ed G a,b i A4 ) 5 


* 
rhus 4 must be the time-symmetric Green function, and can 


cake i. 
be written: G = 4 re C r3 adv where Gret | 


are the retarded and advanced Green functions. 


and. E 


By requiring that the action be stationary under variations 


of the Wisa Hoyle and Narlikar obtain the field-equations: 
ل ظ‎ fa) b Aw ¥ 1 Q 

E 6 M ( X) (x) ] R 2 Qik K) 

aF b 
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wWiere ml) z (G^(x,a) da . However, as a 


consequence of the particular choice of Green function, the 


contraction of the field-equations is satisfied identicalli;. 


‘here ave thus only 9 equations for the 10 components of qu 


ا س 
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and the system is indeterminate. ) 
(o. 


e 
Hoyle and Narlikar therefore impose fm, = ^, =const., 


aS the tenth equation. By then makin» the 'smooth-fluid' 
(a) lb ! 
approximation, that is by putting ŽE ja بي‎ o m 
0. +6 | i 
they obtain the Einstein field-equations: 


C ( . e 
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Thereis an important difference, however, between thes 


field-equations in the direct-particle interaction tceory 


In the 


and in the usual general theory of relativity. 


general theory of relativity, any metric that satisfies the 


the field-equations is admissible, but in the direct-particlc 


interaction theory only those solutions of the field-equations 


are admissible that satisfy the additional requirement: 
( 3 © 
m (x) - 2m ex) : £ (x, زه‎ cC a. 
و‎ . 1 , | 
à £(G, (x, a) da t PALEY (x, a) da 
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This requirement is highly restrictive; it will be shown 
thet it is not satisfied for the cosmological solutions o: 
the Einstein field-equations, and it appears that it cannot 
be satisfied for any models of the universe that either 
contain an infinite amount of matter or undergo infinite 
expansion. 

The difficulty is similar to that occurring in 
Newtonian theory when it is recognized that the universe 
might be infinite. 


The Newtonian potential 0 obeys the equation: 


Hé: - Kp (p 70), 


where P is the density. 


In an infinite static universe, e would be infinite, since 
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the source always has the same sign. The difficulty wa: 
ved when it was realized that the universe was expanding, since 
in an expanding universe the retarded solution of the above 
equation is finite by a sort of'red-shift' effect. The 
advanced solution will be infinite by a 'blue-shift' effect. 
his is unimportant in Newtonian theory, since one is free 
to choose the solution of the equation and so may ignore the 
infinite advanced solution and take simply the finite 
retarded solution. 

tients in the direct-particle interaction theory the 


fr, - field satisfies the equation: 


L]m + + im 3 iV (N 20), 


where A is the density of world-lines of particles. às in 

the Newtonian case, one may expect that the effect of the 
expansion of the universe will be to make the retarded solution 
finite and the advanced solution infinite. However, one is 

now not free to choose the finite retarded solution, for the 
equation is derived from a direct-particle interaction action- 
principle symmetric between pairs of particles, and one must 
choose for fv. half the sum of the retarded and advanced 
Solutions. We would expect tois to be infinite, and this is 


Shown to be so in the next section. 


4. The Cosmological solutions 


The Robertson-Walker cosmological metrics have the 


form . | 4 . ap. T 
$ 2 a PE. Ud 0 | 
io = ct .R* 4 de^. + r ct 0 f SH 

| Ll = AF li. 
Since tney are conformally flat, one can choose coordinates 


in which they become 3 | A 
ds? . JY [as dp. j eles f o siga d 2), 
4 
2 eL oe 
pp 2 1 dx“ 


where Jah io the flat-space metric tensor and 
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For example, for the &Einstein-de Sitter universe 
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For the steady-state (Ge Sitter) universe 
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The Green function C" Co. b) 
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From this it follows that 
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obeys the equation 
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Jg $ * (a, b) 
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This is simply the flat-space vr^en function equation, and 


nence 


q (5,0; t, p) - : 


The '^.-fielà is given by 


men) =(G"W/—g dac 


Por universes without creation (e.g. 


V * 4 FL. 3 


universe), 


Sez 5.) 
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the Binstein-de Sitter 


fe = const. For 
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universes vith creation (steady state) IV fL. fi. 3 const. 
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where the integration is over the future light cone. TALS 

will normally be infinite in an expanding universe, e.g. in 

the isinstein-de sitter universe. 
29 nO 
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In the steady-state universe 
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where the integration is over the past light cone. This will 


normally be finite, e.g. in the Ezinstein-de sitter universe 
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wüile in the steady-state universe 
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Thus it can be seen that the solution M = const. of the 


equation 


| 
im bbe W 


is not, in a cosmological metric, the half-advanced plus 
half-retarded solution since this would be infinite. In fact, 
in the case of the Hinstein-de sitter and steady-state metrics, 


it is the pure retarded solution. 


4, The 'C'-Pield 


Hoyle and Narlikar derive their direct-particle 


interaction theory of the 'C'-field from the action 


A. ££ [{ la, Dui da old! 


where the suffixes Q, 7 refer to differentiation of 


AN 
G (a,b) on the world-lines of (^, b نوم‎ 7 
N 


6 is a Green function obeying the equation 


né x) = £O, X) 
= 
we define the 'C'-field by 


C(x) = Z [G 9); da, 


and the matter-current Jg by | 
7 (y) 2| 5 | 
Then C (x ) : [d^ (x, 4) J Cy); 7 /-9 dx 


IC - ل‎ J K 


ve thus see that the sources of the 'C'-fiela are the places 
where matter is created or destroyed. 


iS in the case of the 'wv'-field, the Green function 


IX x 


must be time-symmetric, that is 


Q (a,b) 3 L ¢ (o, b+ à Qu, CS 4) 


| tet 


Hoyle and Narlikar claim thet if the action of the 


'G'-fielá is included along with the action of the '™ '-field, 
a universe will be obtained that approximates to the steady- 
state universe on a large scale although there may be local 
irregularities. In this universe, the value of C will be 
finite and its gradient time-like and of unit magnitude. 


Given this universe, we may check it for consistency oy ظ‎ 


claculating the advanced and retarded 'C'-fields and finding 
if their sum is finite. We shall not do this directly ^ut 
will show that the advanced field is infinite while the retar- 
ded field is finite. 


Consider a region in space-time bounded by a three- 


dimensional space-like hypersurface D at the present time, 


of » . 


By Gauss's theorem 


AC /-g dx" = ald 


4 n 
£rVv 
1 ١ 
> name J- a x 
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Let the advanced field produced by sources within V be C 5 
2 j 4 Th i 
‘Then. C and will be zero on PA , and hence 


dn 


TK 
But d is the rate of creation of matters N (const.) in 
JK 


the steady-state universe, and hence 


Í oC is 2 n. V -; 
7 dn 


59 ; TJ s " i T i 5 
AS the point F is taken further into the future, the volume 
of the egion V tends to infinity. However, the area of the 
D 
hypersurface tends to a finite limit owing to horizon 
effects. Therefore the gradient (e must be infinite. 
On 
A Similar calculation shows the gradient of the retarded 
field to be finite. ‘their sums cannot therefore sive the 
field of unit gradient required by the Hoyle-Narlikar theory. 


It is worth noting that this result was obtained 


without assumptions of a smooth distribution of matter or 


of conformal: flatness. 


5, Conclusion 


با 


the Einstein theory of 


aJ 
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It is one of the weaknesses Oo: 
relativity that although it furnishes field equations it does 
not provide boundary conditions for them. Thus it does not 
prive a unique model for the universe but allows a whole series 
of models. Clearly a theory thet provided boundary conaitions 
and taus restricted the possible solutions would be very 
attractive. The Hoyle-Narlikar theory does just thst(tne 
requirement thet M = 2% ان‎ 5 1 18 
equivalent to a boundary condition). Unfortunately, as we 
have Seen above, this condition excludes those 700615 that 
seem to correspond to the actual universe, namely the 
Robertson-walker models. 

The calculations given above have considered the universe 
as being filled with a uniform distribution of matter. “his 
is legitimate if we are able to make the 'smootn-fiuid' 
approximation to obtain tne tiinstein equations. Alternatively 
if this approximation is invalid, it cannot be said that the 
theory yields the Einstein equations. 

It might possibly be that local irregularities could make 

WE d finite, but this has certainly not been demonstrated 
and seems unlikely in view of the fact that, in tne Hoyle- 


ct | £2 7 
U -field, 


Narlikar direct-particle interaction theory of their 


which is derived from a very similar action-principle, it can 


be shown without assuming a smooth distribution that the 
advanced 'C' field will be infinite in an expanding universe 
| witn creation. 


The reason that it is possible to formulate a direct- 


particle interaction theory of electrodynamics that does not 
encounter this difficulty of having the advanced solution 
infinite is that in electrodynamics there are equal numbers 
of sources of positive and negative sign. ‘Their fields can 
cancel each other out and the total field can be zero apart 
from local irregularities. This sugzest that & possible way 
to save the Hoyle-Narlikar theory would be to allow masses of 
both positive. and negative sign. The action would be 

A -= ££. 9, foe b) da clb (def b = 5] 


o. Lb 


- 


where Jas (h are gravitational charges analogous to 
electric charres. Particles of positive T in & positive 
'm '-field and particles of negative 4“ in a negative ';^ '- 


field would have the normal gravitational properties, that is, 


toey would have positive gravitational and inertial masses. 


MM Das 000 | | | | 0 nl 


ee 


A particle of negative 1 in a positive 'm'-field would 
still follow a geodesic. Therefore it would be attracted by 
a particle of positive CG . its own gravitational effect 
however would be to repel all other particles. ‘hus it would 
nave the properties of the negative mass described by Bondi 9) 
that is, negative gravitational mass and negative inertial 
Mass. 

Since there does not seem to be any matter having 
tnese oroperties in our region of space ( where m 4A const. > © ) 
there must clearly be separation on a very large scale. It 
would not be possible to identify particles of negative pA 
with antimatter, since it is known that antimatter has positive 
inertial mass. Uwever, the introduction of negative masses 


would probably raise more difficulties than it wuld solve. 
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CHAPTER 2 


PERTURBAT IONS 


1. Introduction 

Perturbations of a spatially isotropic and homogeneous expanding 
universe have been investigated in a Newtonian approximation by 
B... UÁÀ and relativistically by Lifshitz 2), Liftshitz and 
BE atnikov'>) and Irvine M, Their method was to consider small 
variations of the metric tensor, This has the disadvantage that the 
metric tensor is not a physically significant quantity, since one 
cannot directly measure it, but only its second derivatives. It is 
thus not obvious what the physical interpretation of a given 
perturbation of the metric is. Indeed it need have no physical 
significance at all, but merely correspond to & coordinate trans- 


formation, Instead it seems preferable to deal in terms of 


perturbations of the physically significant quantity, the curvature. 


2, Notation 

Space-time is represented as a four-dimensional Riemannian space 
With metric tensor Sap of signature +2.  Covariant differentiation 
in this space is indicated by a semi-colon. Square brackets around 
indices indicate antisymmetrisation and round brackets symmetrisation. 


The conventions for the Riemann and Ricci tensors are:- 


Mec; f bc) - EP ab Vy 0 
R «b + Ra gg 


B. ممع‎ is the alternating tensor, 


Units are such that k the gravitational constant and c, the speed of 


light are one. 


3, The Field Equations 


We assume the Einstein equations: 
Rab = = Gab R 2 ab 


where T nb is the energy momentum tensor of matter, We will assume 


that the matter consists of a perfect fluid. Then, 


Fab = ۸ Ua Uy 4 h hob 


where 0 ls the velecity of the fluid, 8 H* a «4 
w is the density . 
h is the pressure 
Rob = C ou * Ua M is the projection operator 


into the hyperplane orthogonal to U, 


hab u” 2 ©, 


We decompose the gradient of the velocity vector U, as 


where TF = Uab “ابا‎ is the acceleration, 

D = u,'* is the expansion, 

O,& = (e;d) he hí -$ hap is the shear, 

Gap = reid] كط‎ hy is the rotation of the 
flow lines .وا‎ We define the rotation vector W, as 

Wa * € Jabed o»** uu? 


We may decompose the Riemann tensor B. ond into the Ricci tensor 


R., and the e 8 | 
ab und the Veyl tensor G abod’ 


R abed = Cabed - Sopa Rc] b = Supe Raga ~ P Sals Je) b 
C abed 2 C (o. b] رمع‎ 


) 


hen SOA 4 of be 


C abcd 


is that part of the curvature that is not determined locally by 


the matter. lt may thus be taken as representing the free gravit- 


ational field (Jordan, Ehlers and romat, We may decompose it 


into its "electric" and "magnetic" components. 


= = 5 
Es 2d C <b paq uF U 2 


c d - L6 3 d 
C e oi 5 8 Ure Ep ud - 4 pa ET à 


d cdr 
-2 Nabea U H^l*u*l. om acd 7 Here 


|- = € | - 
Zak (eb) د‎ Ha = Hab , 


E and Ho each have five independent components. 


ab p 
We regard the Bianchi identities, 


et way eo’) 59 


as field equations for the free gravitational field. 


Then 


sal 
C abed^" = + Retos] Sepp Ra] (Kundt 


vi 
b) 2» 


and Trümper, 


*6)), 


Using the decompositions given above, we may write these in a form 


analogous to the Maxwell equations. 


de 


cd c 
b. E bcd h + 5 Haso” - 7) طن‎ cd uo eH - ha” usb 3 (1) 


c a ee 
h. ^H be: d "n - 3 Fab w? "- YJabcol u^g e = a (Ki h) Wa , (2) 


—— 


| | , 4 d;e c 
IX E. t ho Yede ut Hye + EO - E (ey. 
-E (a Mb). ^ Yocde Y bpqr uo E: E*' 


d 6 a: í 
+2 H'a peas عنا علا‎ = - (Mth) Tab ; (3) 


5 
E na- higa" ule W E, ap Hag = H ea W bye 


H (a Te - Noede مومه[‎ U uP ON H 77 (4) 
FE Hells us ue م‎ r 
where l indicates projection by سمط‎ orthogonal to U.. 
io. f. ال‎ 
The contracted Bianchi identities give, 
/ a nee i 
(Ran ¥ 4abR) = a Pm ) 
p + 06ج مر)‎ = O , (5) 
b 
(«+ h) lig + Rp H's ماع‎ : (6) 


The definitien of the Riemann tensor is, 


U aib] = 2 Rapbe UP 
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Using the decompositions as above we may obtain what may be regarded 


as "equations of motion", 


7 2 MERE" 
B -ac'-20'- $0 + ua - v(pest) , 


o : 1 3 q 
10 لف‎ ob = ti $ Wale 9 ف‎ 2 Oeta (Day + لما‎ £p;a] hc h b 
| E 
| 1 
| £L O ab " E ab 3 GO ud % ل ب‎ 0% - + Og, 


| م‎ s e : 
udi. 1 ak iG s atg eae T Ug Ab 


+ كط ورمعن‎ hl 


2 Ly 2 as OD 
where 2 69" ب‎ Web GO 8 20 = Oab O 


We also obtain what may be regarded as equations of constraint. 


guh = 2 (eyes o Pee) كط‎ ah Qoo + Sat) 


* & i a. 


W) م‎ = 2 Wau ) 


d:‏ - كن 
Hu = hee u* [co ,7* + 903 ie | à‏ 


(10) 


(11) 


(12) 


We consider perturbations of a universe that in the undisturbed state 


in conformally flat, that is 


U سم‎ s o 8 


By equations (1) - (3), this implies, 


Gab = Wab= 0 


i b 
ha? رهز‎ = oc 6.4, ha 


If we assume an equation of state of the form, f= (PM) , 


{nen by (6), (10), hyo ha = Q 2 in 4 


This implies that the universe is spatially homogeneous and isotropic 
since there is no direction defined in the 4-space orthogonal to Use 

In this universe we consider small perturbations of the motion 
of the fluid and of the Weyl tensor. We neglect products of small 
quantities and perform derivatives with respect to the undisturbed 
metric. Since all the quantities we are interested in with the 
exception of the scalars, Es Te, O have unperturbed value zero, we 
avoid perturbations that merely represent coordinate transformation 
and have no physical significance. 


To the first order the equations (1) - (4) and (7) - (9) are 


Ja (13) 
Hay” = (pth) Wa ; (11) 
E au + ESO £^ eoe U HUS. -z (pth) Tab» (15) 
Hab *HasB د‎ h^. Yu ade ut جح‎ S we. (46) 
6 ب‎ FO + و ي5#+سن) د ی‎ (17) 
Quac TE ab 0 + Agea 5658 , ban 
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e] nA "i Pt 
O ob " E b i = 6,, 0 ~ 9 rts Ue 5 T رمع لما‎ q) n. h (19) 


from these we see that perturbations of rotation or of E b or Hop do 


not produce perturbations of the expansion or the density. Nor do 


perturbations of E b and Hoy produce rotational perturbations. 


h, The Undisturbed Metric 
Since in the unperturbed state the rotation and acceleration 


are Zero, 8 must be hypersurface orthogonal. 


where T measures the proper time along the world lines. As the 
surfaces T = constant are homogeneous and isotropic they must be 
4esurfaces of constant curvature. Therefore the metric can be 


written, " 
ds" z-dct*a.4 $t dy” 


re ç? - Qv) ? 
dy is the line element of a space of 


zero or unit positive or negative curvature, 


We define t by, 


dt. ل‎ 
At i 
2 
then إن‎ s* 3 «At i + dy”) 1 
In this metric, Ua = ( - 9 ME &, Os, o) 


, o = 3560 . BY 


9 9 


(prime denotes differentiation with respect to t) 


Then, by (5), (7) 


b=- lh) i و‎ 


(20) 
z C) l ( J 
هد ج‎ wee i A T5 r) 5 2۹ 
c / | (21) 
If we know the relation between u and fi, we may determine (l i 
We will consider the two extreme cases, f^ - O (dust) and h- A 
(radiation). Any physical situation should lie between these. 
Por f= 9 
By (20) - M Nk = const 
9 F : a. 4 = o 
* i 
3 f "m 
due 2h. < = 6 
i M y :£* 
: 5 ge. whee 
mA 6 E ; E - const. 
(a) For E O, 
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(c) For E O, E 
$2 = zE (! cos J OEC t) , 2E NER , 


E represents the energy (kinetic + potential) per unit mass. 


If it is non-negative the universe will expand indefinitely, other- 


Wise it will eventually contract again. 


By ine 


Gauss Codazzi equations *R,the curvature of the 


hypersurface T = const. is bom of 4 na ق‎ ) 
IN e 4 2 v J 
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5. Rotational. Perturbationg 
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Thus rotation dies away as the universe expands. This is in fact a 


statement of the conservation of angular momentum in an expanding 


universe. 


6. Perturbations of Density 


For A. = Ó we have the equations, 


i 


Mox TAE 


an di ا‎ Ta 


1 
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These involve no spatial derivatives. Thus the behaviour of one 
region is unaffected by the behaviour of another. Perturbations 
will consist in some regions having slightly higher or lower values 
of E than the average. If the universe as wh*te has a value of E 
greater than zero, a small perturbation will till have E greater 


than gero and will continue to expand. It will not contract to 


form a galaxy. If the universe has a value of E less than zero, a 


ت 


small perturbation can contract. However it will only begin 
contracting at a time 8T earlier than the whole universe begins 


contracting, where 


C ^T^ aM. sa 
50 |. DG 
و‎ — " - 
VE 5 m 
T" bs 


© o is the time at which the whole universe begins contracting. 
There is only any real instability when E = O. This case is of 
measure zero relative to all the possible values E can have. 
However this cannot really be used as an arguement to dismiss it 


as there might be some reason why the universe should have it = O, 


For a region with energy -8E و‎ in a universe with E = O 


For E = O, Ms tt 


Thus the perturbation grows only as T . This is not fast enough 
to produce galaxies from statistical fluctuations even if these 
could occur. However, since an evolutionary universe has a particle 
horizon Hindle, Panos i?) different parts do not communicate 
in the early stages. This makes it even more difficult for 


Statistical fluctuations to occur over a reg:'n until light hed time 


to cross the region, 
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as before, a perturbation cannot contract unless it has a negative 
value of E. The action of the pressure forces make it still more 


difficult for it to contract. Eliminating O, 
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to our approximation. 
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These eigenfunctions will be hyperspherical and pseudohyperspherical 
harmonics in cases (c) and (a) respectively anc plane waves in case 


(b), In ease (c) n will take only discrete ` (lues but in (a) and 


(b) it will take all positive values. 
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where Ho is the undisturbed density. 
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These perturbations grow for as long as light has not had time to 
travel a significant distance compared to the scale of the perturbation 
E fl ). Until that time pressure forces cannot act to even out 


perturbations. 
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"Ue obtain sound waves whose amplitude decreases with time. These 
results confirm those obtained by Lifshitz and Khalatnikov ^). 
From the forgoing we see that galaxies cannot form as the result 


of the growth of small perturbations. We may >xpect that other non- 


gravitational force: will have an effect smaller than pressure equal 


to one third of the density and so will not cause relative perturbations 


to grow faster than T . To account for galaxies in an evolutionary 
universe we must assume tnere were finite, non-statistical, initial 


inhomogeneities. 


| 7. The Steady-state Universe 


To obtain the steady-state universe we must add extra terms to 
(10) use, 


the energy-momentum tensor, Hoyle and Narlikar 
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There is a difficulty here, if we require that the "C" field 


ghould not produce acceleration or, in other words, that the matter 


ereated should have the same velocity as the matter already in 


oxistence,We must thën have 


he. Ck s5 0 


However since C is a scalar, this implies that the rotation of the 
medium is zero. On the other hand if (23) does not hold, the equations 
ape indeterminate (c.f. Raychaudhuri and Bernop4eg! 123, In order to 
have a determinate set of equations we will adopt (25) but drop the 
requirement that C, is the gradient of a scalar. The condition (23) 
is not very satisfactory but it is difficult to think of one more 
Satisfactory. Hoyle and Narlikar V?) seek to avoid this difficulty 

by taking a perticle rather than a fluid picture. However this has a 
serious drawback since it leads to infinite fields adds 243, 
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Thus, small perturbations of density die away. Moreover equation (18) 


still holds, and therefore rotational perturbations also die away. 
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These results confirm those obtained by Hoyle and Narlikar We 
see therefore that galaxies cannot be formed in the steady-state 
universe by the growth of small perturbations, However this does not 
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exclude the possibility that there might by a self-perpetuating 


system of finite perturbations which could produce galaxies, 


(15) (6), 


(Sciama , Roxburgh and Saffman 
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We now consider perturbations of the Weyl tensor that do not 


arise from rotational or density perturbations, that is, 


Multiplying (15) by US Ve, and (16) by hla Urs; 


we obtain, after a lot of reduction, 
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In empty space with a non-expanding congruence U this reduces to 


the usual form of the linearised theory, 
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The second term in (2h) is the Laplacian in the hypersurface 
0] = constant, acting on Ey We will write Ep as a sum of 


eigensfunetions of this operator. 
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80 the gravitational field E , decreases as 91 7! and the "energy" 
(Bap ga Hp He? ) as p; » We might expect this as the 


Bianchi identities may be written, to the linear approximation, 
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Therefore if the interaction with the matter could be neglected 
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ica would be proportional to © and Bh f Lp to (1 ‘ 
In the steady-state universe when u and © have reached their 


equilibrium values, R t= 8 + A) ded 
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mous the interaction of the "C" field with gravitational radiation is 


equal and opposite to that of the matter. There is then no net 


* 1 A ^ 
interaction, and E, and H decrease as [OQ ' 
ab ab 
1 m e ab TT ab z. 0 0 
The "energy" 2 (E pË HH ) depends on seeond derivatives 
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of the metrie. It is therefore proportional to the frequeneéy squared 
times the energy as measured by the energy momentum pseudo-tensor, in 
a local co-moving Cartesian coordinate system which depends only on 
first derivatives, Since the frequency will be inversely proportional 
to (2 و‎ the energy measured by ihe pseudo-tensor will be proportional 


to C H as for other rest mass zero fields, 


9. Absorption of Gravitational Waves 

As we have seen, gravitational waves are not absorbed by a 
perfect fluid. Suppose however there is a small amount of viscosity. 
We may represent this by the addition of a term A Cab to the 
energy-momentum tensor, where is the coefficient of viscosity 
ers, 17). 


Since Tab’ =o 


we have 


qquations (15) (16) become 
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The extra terms on the right of equations (27), (28) are similar to 

conduction terms in Maxwell's equations and will cause the wave to 
VE 

decrease by a factor e Ze . Neglecting expansion for the moment,. 


suppose we have a wave of the forn, 


This will be absorbed in a characteristic time “> independent of 
frequency. By (25) the rate of gain of rest mass energy of the 
matter will be 2^4 ^ which by (19) will be 2A E v. Thus the 
available energy in the wave is  *4 QC Y. This confirms that the 
density of available energy of gravitational radiation will decrease 
ag [7 s: in an expanding universe. From this we see that 
gravitational radiation behaves in much the same way as other 
radiation fields. In the early stages of an evolutionary universe 
When the temperature was very high we might expect an equilibrium to 


bo set up between black-body electromagnetic ; diation and black-body 


gravitational radition. Since they both have two polarisations their 


energy densities should be equal. As the universe expanded they woulc 


poth cool adiabatically at the same rate. As we know the temperature 
of black-body extragalactic electromagnetic radiation is less than 


BOX , the temperature of the black-body gravitational radiation must 


be also less than this which would be absolutely undetectable. Now 
the energy Cf gravitational 2503 tion does not contribute to the 


ordinary energy momentum tensor Nevertheless it will have an 
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active gravitational effect. By the expansion equation, 
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G 
radiation has an active attractive gravitational effect. It is 
interesting that this seems to be just half that of electromagnetic 
radiation. 


and Hoyle and Nariikse O, 


(18) 


It has been suggested by Hogarth 
that there may be a connection between the absorption of radiation 
and the Arrow of Time. Thus in universes like the steady-state, in 
Which all electromagnetic radiation emitted is :ventually absorbed by 


Other matter, the Absorber theory would predic retarded solutions of 


whei u is the gravitational "energy" density. Thus gravitational 


the Maxwell equations while in evolutionary universes in which 


electromagnetic radiation is not completely absorbed it would predict 
advanced solutions, Similarly, if one accepted this theory, one would 
expect retarded solutions of the Einstein equations if and only if 3 
gravitational radiation emitted is eventually absorbed by other matter, 
Clearly this is so for the steady-state universe since ^ will be 
constant, In evolutionary universes ^ will be a function of time. 

We will obtain complete absorption if \ Ade diverges. Now fr? a gas, 
TX where T is the temperature. For a monatomic gas, TEN’, 
therefore the integral will diverge (just). However the expression 

used for viscosity assumed that the mean free path of the atoms was 
small compared to the scale of the disturbance. Since the mean free 
path عه‎ ul (27? and the wavelength « (1^! , the mean free path will 
eventually be greater than the wavelength and so the effective viscosity 
will decrease more rapidly than (2 » Thus there will not be comolcte 
absorption and the theory would not predict retarded solutions. 

However this is slightly academic since gravitational radiation has nol 


yet been detected, let alone investigated to see whether it correspmds 


to a retarded or advanced solution, 
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CHAPTER 2 


Gravitational radiation In An 


Expanding Universe 
Gravitational radiation in empty asymptotically flat 


Space has been examined by means of asymptotic expansions 
[4] 


by a number of authors. "They find that the different 
components of the out;oing radiation ficld "peel off", that 
is, they go as different powers of the affine radial distance. 
If one wishes to investigate how this behaviour is modified 
by the presence of matter, one is faced with a difficulty 
that does not arise in the case of, say, electromagnetic 
radiation in matter. For this one can consider the radiation 
travelling through an infinite uniform medium that is static 
apart from tne disturbance created by the radiation. In the 
case of gravitational radistion this is not possible. tor, 
if the mecium were initially static, its own self sravitation 
would cause it to contract in on itself and it would cease to 
be static. Hence one is forced to investigate gravitational 
raalation in matter that is either contracting or expanding. 
48 in Chapter 2, we identify the Weyl or conforma 
tensor NE with the free gravitational field and the 


- , () . 1 
Ricci-tensor hal with the contribution of the matter to the 


curvature. Instead of considering gravitational radiation in 


tically flat space, that is, space that approaches 


flat space at large radial distances, we consider it in 
asymptotically conformally flat space. As it is only 


conformally flat, the Ricci-tensor and the density of matter 


| need not be zero. 

To avoi&à essentially non-gravitational phenomena such 
as sound waves, we will consiaer gravitationalraciasvion 
travelling through dust. It was shown in Chapter 2 that a 
conformally flat universe filled with dust must have one of 
the metrics: 
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"ype(a) represents a universe in which the matter 
expands from the initial singularity with insufficient energy 
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to reach infinity and so falls back again to another 
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Singularity. It is therefore unsuitable for a discussion of 
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gravitational radiation by a method of &symototic expansions 


KES 


since one cannot set an infinite distance from thi@ source. 

Type (b) is the Zinstein-De Sitter universe in which the 
matter nas just sufficient energy. to reach infinity. It is 
thus a Special case. D. Norman (2) has investigated the 
"peeling off" behaviour in this case using Penrose's conformal 
technique (6). He was however forced to make certain assumpt- 
10215 about the movement of the matter which will be shown to 
be false. Moreover, he was misled by the special nature of 


the Hinstein-De Sitter universe in which affine and Luminosity 


HP. 
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stances differ. Another reason for not considering radiation 
in the isinstein-De Sitter universe is that it is unstable. 
Phe passage of a gravitational wave will cause it to contract 
azain eventually and develop a singularity. 

4e will therefore consider radiation in a universo of 
type (c) which corresponds to the general case where the 
matter is expanding with more than enough energy to avoid 


contracting again. 


2. ‘the Newman-Penrose Formalism 


ve employ the notation of Newman and Pantoss, 4 A 


I^ m^ = 
tetrad of null vectors ,U, IT , Ifi m^ is introduced 
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azain eventually and develop a singularity. 

4e will therefore consider radiation in a universe of 
type (c) which corresponds to the general case where the 
matter is expanding with more than enough energy to avoid 


contracting again. 
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we label these vectors with & tetrad index 
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tetrad indices are raised ana lowered with the metric 
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Ricci rotation coefficients are defined by: 
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In fact it is more convenient to work in terms of twelve 
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complex combinations of rotation coefficients defined as 


follows: 
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2, Coordinates 


Like Newman and Penrose, we introduce a null coordinate 
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geodesic and irrotational. This implies 
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we take I^ m” -m to be parallelly transported 
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along JF . This gives 
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AS a second coordinate we take an affine parameter (= x J 


long the geodesics 4 
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In these coordinates 


The Field Ecuations 
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We may calculate the Ricci and Weyl tensor components from 
the relations elc 
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Using the combinations of rotation coefficients already 


defined and with Ke TT = €.20 we have 
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Expressin; the rotation coefficients in terms of the metric, 


we have: 


yr das d s e3) 
DO د‎ gos eo-(G«g) QOD 
DX DE) 1 E G68) 
DU > Ò= ïo- lysg) Q9 
KA Ge poa) (G4 
pei .نمق‎ GE GPE Gs 
ga -So = (B -a)or G-G)5sQe-F) Éu) 
$U- Ao: (py - p)77 \a-% G59 


As in Chapter 2 we use the Bianchi identities as 


field equations for the Weyl tensor. In tne Newman-renrose 
formalism they may be written: 


(I am indebted to R. G. McLenaghan for these) 
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4, The Undisturbed Metric 
The undisturbed metric may be written 

42 1 | "4. ١ 2 2. i bz) 
ds? : 3L Cae ad gt sinh e ) Ao rSin d | 

$4. 5 4 (cosh & =1) 
put P عد‎ o -0 

| | lef, IRE. Bt كن‎ 

NE do? 1 m d dix? 3 JA ue - siah (t- u Gy sin Clg )| 
(4 0) 


Wù is a null coordinate 
To calculate وم‎ the affine parameter, we note thet C 
is an affine parameter for the metric within the square 
brackets. Therefore = 121 td 6s («OG p) (4.2 ) 
will be an effine parameter for (“1) 
]j is constant along the null geodesic. Normally it 
would be taken so that [f= O when E i . However, 


in our case it will be more convenient to make it zero and 
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This means that surfaces of constant f ape surfaces of 


define i as 


constant C . This may seem rather odd, but it should 
be pointed out that the choice of X will not affect the 


asymptotic dependence of quantities. That is, if 
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Then 


It proves easier to perform tie calculations with this 


choice of f but all results could be transformed back 


to a more normal coordinate system. 
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The matter in the universe is assumed to be dust so its 
energy tensor may be written 
ete V (4.5) 
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For the undisturbed case, from Chapter 2 
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Now 
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Therefore if we try to expand (۸ as a Series in powcr of S 
the result will be very messy and will involve terms of 

the form log S * 

*Tt should be pointed out that the expansions used will 
only be assumed to be valid asymptotically. They will not 
be assumed to converge at finite distances nor will the 


quantities concerned be assumed analytic. (see A. Erdelyi: 


Asymptotic Expansions - Dover 
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This Goes not invalidate it as an asymptotic expansion but 


it makes it tedious to handle. For convenience therefore, 


we will perform the expansions in terms of err ts which 
will be defined in general as the same function of j^ as 
it is in the undisturbed case. That 1s i 


where ~ = D> [ + se^ 2€ - Decr EC 


unen 


| 
: 


! (d م‎ a SA oR? 26 
|] * ure aj 1 y wv ف‎ Gl 

— | 2ES a 55S 

$& 23? QS ey) $2 - 
For the third and fourth coordinates it is more convenient to 
use stereographic coordinates than spherical polars. 

Since the matter is dust its energy-momentum tensor 

and hence the Ricci-tensor have only four independent 
components. We will take these as / D , Qs , 

. / 100 í 
(Since Do, is complex it represents two components) 
In terms of these the other components of the Ricci-tensor 
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For the undisturbed universe with the coordinate system 


given: A = po 9 A | 
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Using these values and the fact that in the undisturbed 
universe all the ye are zero, we may integrate equations 
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5. Boundary Conditions 


"e wish to consider radiation in a universe that 


asymptotically approaches the undisturbed universe given 


above. on and /\ will then have the values 
ziven above plus terms of smaller order. To determine this 
order and the order of (D. and y5 yg there- are two 
ways in which we may proceed. we may take the smallest orders 
chat will permit radiation, that is V, T Cr U . 
Larger order terms than these in D, 9 A and. 

cp turn out to have their UL derivatives 
BS 2a ant only on themselves and not on the r^ coefficient 
of V, , the radiation field. ‘hey are thus disturbances 


not produced by the radiation field and will not 26 


Alternatively we may proceed by a method of successive 


approximations. We take the undisturbed values of the spin 


coefficients and use them to solve the Bianchi Identities as 
field equations for the conformal tensor using the flat space 
boundary condition that LE ole?) . Then substituting 
these y" in equations (3.jO- 4F ) calculate the disturb- 
ances inauced in the spin coefficients and substituting these 
back in the Bianchi Identities, calculate the disturbances in 
the ap" . ij'urther iteration does not affect the orders 
of the disturbances. 
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it will be shown that if these boundary conditions 
hold on one hypersurface (ù = const.) they will hold on 
succeeding hypersurfaces and that these conditions are the 


most severe to permit radiation. 


Integration 
AS Newman and Penrose, we begin by integrating the 
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Tous the tw derivative of b. ; depenas only on itself 

and not on the radiation field. lt therefore represents a 
type of disturbance unconnected with radiation. If it is 
zero on one hypersurface, it will remain zero. In this case 
it is possible to continue the exoansions of ail auantities 


in negative powers of JA without any log terms appearing. 
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The asymptotic group is the group of coordinate transformations 
that leave the form of the metric and of the boundary conditions 
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che asymptotic group is isomorphic to the conformal group in 
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| to the group of motions of a -space of constant negative 
curvature which is the group of the unperturbed Robertson- 
;alker space. Thus the asymptotic group is the same as the 

| group of the undisturbed space. It is not enlarged by the 


resence of radiation. This is interesting because in tho 


um 


i. 


case of gravitational radiation in empty, asymptotically flat 
pace, it turns out that the asymptotic group contains not 
only the 10 dimensional inaomozeneous Lorentz group, the group 


of motions of flat space, but also infinite dimensional | 


"Supertrauslations". It has been suggested that these 
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lf the Linstein equations without cosmological constant 
are satisfied, a Robertson-Walker model can 'bounce' or avoid 
a singularity only if the pressure is less than minus one- 
third the density. This is clearly not a property possessed 
by normal matter though it might be possessed by a field of 
negative energy density like the 'C' field. However there is 
a grave quantum-mechanical difficulty associated with the 
existence of negative energy density, for there would be 
nothin, to prevent the creation, in a given volume of space- 
time, of an infinite number of quanta of the negative enersy 
field and a corresponding infinity of particles of positive 
energy. lf we therefore exclude such fields, all Robertson- 
Walker models must be of the 'big-bang' type, that is they 
nave a Singularity in the past and maybe one in the future 
as well. lt has been suggested k that tne occurrence of 
these singularities is a consequence of the high degree oi 
symmetry of the Robertson-Walker models which restricts 6 
expansion and contraction so that they are purely radial and 


that more realistic models with fewer or no exact symmetries 


would not have a singularity. This chapter will be devoted 


to an examination of this question and it will be snown that 


provided certain physically reasonable conditions hold, any 
model must have a singularity, that is, it cannot be a 
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als 
Therefore d (Fl) is finite 
oL Mn 
3 
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Hence & is infinite where and only where clol Mi 2 0 
Thus the two definitions of conjugate points are equivalent. 
This also shows that singular points of congruences. are 


points where neighbouring geodesics intersect. 
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lor null geodesic congruences with parallelly transported 


17 


tangent vector ^de) we may define the convergence e Qs 
in Chapter j. this obeys ۴ 
Bat = are + ak, See Sy (s) 
we define a singular point of a null geodesic congruence as 
one where d is infinite. 
whe condition that tne pressure is greater than minus 
one-third tule density may be stated more generally as 
condition (a). 
(aj pow. E ^X i for any 
observer with 4-velocity 7 where -= neun 
is the energy density in the rest-frame of the observer and 
why ١ 9 ls the rest-mass density. 
Condition (a) will be satisfied by a perfect fluid with | 


. "Ic 
density حمر‎ 2 © and pressure P2 Kd lt implies kV |/ 20 for 


any time-like or null vector y^ e ‘therefore by equations 

(1) and (5) any time-like or null irrotational geodesic 

congruence must have a singular point on each geodesic within 

a finite affine distance. Obviously if the flow-lines form 

on irrotational geodesic congruence, there will be a physical 

Singularity at the singular points of the congruence where 

the density and hence the curvature are infinite. ‘this will 
| 


be the case if the universe is filled with non-rotating dust 


However, if the flow-lines are not geodesic (ie. non-vanishing 


pressure gradient) or are rotating, equation (1) cannot be 
applied directly. 


atially ilomozeneous Anisotropic Universe 


io 


The Robertson-wW/alker models are spatially homogeneous 


and isotropic, that 1s, they have a six parameter group of | 
motions transitive on a sp&ce-like surface. lf we reduce the 
symmetry ùy considering models thét are spatially homogeneous 


but anisotropic (that!is, they have a three parameter group 
t , J I, c [ 


the matter flow may have rotation, acceleration and shear. 
Thus there would seem to be the possibility of non-singular 
models. dL. 8 epley' has investigated one particular 
homogeneous model containing rotating dust and has shown that 
there is always a singularity. Here a general result will be 
proved. 

there must be a Singularity in every model which satisfies 
condition (a) and, 
(b) there exists a Gy of motions on the space or on 


universal covering space  , P S which is transitive on at 


* See section 5 


— — — — 


YO 


least one space-like surface but space-time is not stationary, 
(c) the energy-momentum tensor is that of a perfect fluid, 

^ i U, لا‎ - Paap . U^ is the tangent ic 

. lab س‎ (u row P) C b Gel o UL LS CHC UA CNH U U U 
flow-lines and is uniquely defined as the time-like eigen- 


vector of the kRicei tensor. 


like surface of transitivity H” of the group.  iherefore Keg 
4 
| oe 


. 1 E . ٠. 1 . LI . 45 
must be in the direction of the unit time-like normal l/ 


to He . 


PROOF 
R, the curvature scalar must be constant on a space- 


RS‏ م 
wnere T 5 f >o‏ 


e lR, ) is an indicator = +1 if dT. is past directed 
pr | 
ب‎ =] if دكا‎ is future directed 
\ : "A 
Lien V. , —1 9 
Lo: b] 
thus Va lS a congruence of geodesic irrotational time-like 


a b 

vectors. By condition (a), - V V 209 i 
therefore the congruence must have a singular point on each 
seodesic ( by equation 1) either in the future or in the past. 
jurther, by the homogeneity, the distance along each geodesic 
م‎ H? + که ع‎ : d. e 1 ; na a 
from to the singular point must be the same. ior eacn geodesic. 
thus if the surfaces of transitivity remain space-like, they 

"E 9 5 ias a i —Áá Sah wil 
must degenerate into, at the most, a 2-surface ) which will 
be uniquely defined . tet M be the subset of the flow-lines 

2 5 

of the matter which intersect ه و‎ Lev B 


3 
empty subset of H intersected by M . wince there is a 


^ 
22 A i 35 ; 
group transitive on H T must be H itself. Thus all the 


D 


7 7 
5 ; x H^ i ks a i "T^. 5 
flow-lines through FH must intersect tne 2-surfaáce C . ‘Thus 


the density will be infinite there anu there will be a 


iysical singularity. Alternatively if the surfaces of 
transitivity do not remain space-like, there must be at 
least one surface which is null - call this S^. At 84. p = QO, 
M y í irk” is zero, we can take any other scalar 
oolynomial in the curvature tensor and its covariant derivatives. 
‘ney cannot all be zero if space-time is not stationery). 3e 
introduce a geodesic irrotational null congruence on 37 with 

Cn i 
tangent vector L where " R. ; ‘Then by equation (5), 

^ / 

there will be a singular point of each null geodesic in 5 
within finite affine distance either inthe future or in the 
past. ‘the 2-surface of these singular points will be uniquely 
defined. The same argument used before shows that the density 


يه 
5 


becomes infinite and there is a physical singularity. in fact 


a 


as 87 is a surface of homoseneity, the whole of S^ will be 


Singular and it is not meaningfulito call it null or to 
distinguish tals case from the case where the surfaces of 
transitivity remain soace-like. 

The conditions (a), (b), (c) may be weakened in two way 


Condition (b) that there is a group of motions throughout 


, , F. — 5 
space-time may be replaced by (b ) and (dû). 


/ 


2 
na ; "a 5 cry 11 ay iii AP es if” TENET CEAEE S^ ر‎ 
ib ) There is a Space-like hypersurface in which there 


it | 


j : J 3 o» S 3 , Or 1 1 
are three independent vector fields A such that 
Mm 


PRU d 
d q bc va K. de = 0 on H 3 ٠. that is, there 
Mm 


15 one homogeneous space section. 
(d) ‘here exist equations of state such that the vaucny 
| 3 . ! 
development of H 1s determinate. 
Then succeeding space-like surfaces of constant R 
are homogeneous and much the same proof can be given that there 


/ 
are no non-singular models satisfying (a), (b), (o), (à). 


intersection of which implies a physical singularity. Obviously, 
however, tais property will be possessed by a much more gener 
class of fluids. For these, we define the flow vector as the 
time-like eigenvector (assumed unique) of the energy-momentuii 
tensor. ‘Then we can replace condition (c) on the nature of 

the matter by the much weaker condition (e). 

(e) If the model is singularity-free, the flow-lines form 
a smooth time-like congruence with no singular points with 


a line tarough each point of space-time. 


Condition (e) will be satisfed automatically if conditions 


(a) and (c) are. 


This proof rests strongly on the assumption of 


homogeneity which is clearly not satisfied by the physical 
universe locally though it may hold on a large enough scale. 
However it would seem to indicate that large scale effec 
like rotation cannot prevent the singularity. 

lt is of interest to examine the nature of the singularity 
in the homogeneous anisotropic models since this is more 
likely to be representative of tne general case than that of 
the isotropic models. It seems that in general the collapse | 
will be in one direction,” that is, the universe will collapse 
down to a e-surface. Near the singularity, the volume will be 
proportional to the time from the singularity irrespective of 
the precise nature of the matter. lt also appears that tne 
nature of the particle horizon is different. ‘There will be 
a particle horizon in every direction except that in which 


the collapse is taking place. 


4, singularities in Inhomozeneous Models 


Lifshitz and Khalatnikov? claim to have proved that a 
general solution of the field equations will not have a 
Singularity. Their method is to contract a solution with a 
singularity which they claim is representative of the 


general solution with a singularity, and then show that it 


has one fewer arbitrary function than a fully general solution. 


Clearly their whole proof rests on whether their solution 


7 


is fully representative and of that they give no proof. 
indeed it would seem that it is not representative since it 
involves collapse in two directions to a 1-surface whereas 
in general one would expect collapse in one direction to 4 
e-surface. In fact their claim nas been proved false by 
Penrose’ for the case of a collansing star using the notion 
of a 'closed trapped surface'. A similar method will be 
used to prove the occurrence of Singularities in 'open ' 


universe models. 
5. 'Open' and 'Closed' Models 


The method used by Penrose to prove the occurrence | 


of a physical singularity depends on the existence of a 
non-compact Cauchy surface. A Cauchy surface will be taken 

to mean a complete, connected spece=like surface that 
intersects every time-like and null line once and once only. 
Not all spaces possess a Cauchy surface: examples of those 
that do not include the plane-wave nelê, the Godel model,” 
and BN.U.T. Space" However none of these have any paysical 
Significance. Indeed it would seem reasonable to demand of 
any physically realistic model that it possess a VUauch 


surface. If the Cauchy surface is compact, the model is 


commonly said to be 'closed'; if non-compact, it is said to 
9 9 


to be 'open'. The surfaces, t = constant, in the Robertson- 


walker solutions for normal matter are examples of Cauchy 
surfaces. If K = -1, they have negative curvature and it is 
frequently stated that they are non-compact. This is not 
necessarily so: there exist possible topologies for which 
they are compact. However, the following statements may be 
made about the topology of the surfaces i = constant. 


If the curvature is negative, X = -1, the universal 
4^ 


WN 
^ 


covering space is non-compact and is diffeomorphic to E 


Any other topology can be obtained by identification of 
points. Thus any other topology will not be simply connected 
and, if compact, must have elements of infinite order.in the 


fundamental group. Further if compact, they can have no 


fo 
; dd 
group of motions. 


If the curvature is zero, K = و0‎ tne universal covering 


sc, ا‎ " eT : | 4 55 
Space is à^. There are eighteen possible topolosies. Í LI 


wT 
PS 


compact they have 2 Gz of motions and Betti numbers, B, = 5, 


l 8 12 
Bo = 2 « 


If the curvature is positive, K = +1, the universal 


i11 


2 
covering space is S^. Thus all topologies are compact. The 
. » 12 
Betti numbers are all zero. 

oince a Singularity in the universal covering space 


implies a singularity in the space covered, Penrose's method 


is applicable not only to spaces that have a non-compact VUauchy 


surface but also to spaces whose universal covering space 


nas a non-compact Cauchy surface. ‘hus it is applicable to 
models which, at the present time, are homogeneous and iso- 
tropic on a large scale with surfaces of approximate homo- 


+ 


geneity which have negative or zero curvature. 


6. The Closed Trapped Surface 
Let p? be à 3-ball of coordinate radius r in a -surface 


H’ (t = const.) in a Robertson-Walker metric with K = © or -1. 


P, 


2 a E1 1 ° ^ m 3 M ا‎ 
Let q be the outward directed unit normal to T^, the boundary 
"amm m 8 pa i ae ee : 
of T7”, in H^ and let V be the past directed unit normal to 


u^. Consider the outgoing family of null geodesics which 


-— — T TC Dus 
intersect T^ orthogonally. at T. و‎ their convergence will 
| 6 | (S^ 9b. t^t) 

f i (Vs. * Garb) 


E 
) a à; i i . 59 : 
where S^ e are unit space-like vectors in H^ orthogonal 


tO a and to each other, 
ji | wie dh. Mcd | 
z CB. -k pyi- kt 


O or -1, by taking r large enough, we may 


! Ks { 2 L 5 Fa. TS 
ai a negative at T`. Therefore, in the language oi Perro se, 


We a 3 " 
T- is a closed trapped surface. 


b 


toerefore ¢ 


if W7O and K 


Another way of seeing tnis is to consider the diagram 


in which the flow-lines are drawn at their proper spatial 


distance from an server. They all meet in ine singularity 


tj 


at t = .نأا‎ 1f the past light cone of the observer is drawn 


on tnis diagram, it initially diverges from his world-line 

( 1 <O ). It reaches a maximum proper radius ( à B0 5 

and then converges again to the singularity ( pou ). The 
2 


intersection of the converging light cone and tne surface H 
gives a closed trapped surface P. If the red-shift of the 
quasi-stellar 209 is cosmological then it will be beyond the 
point e ZO if we are living in a Robertson-walker type 

universe with normal matter. However, the assumptions of 

homogeneity and isotropy in the large seem to hold out to tne 
distance of 5079. ‘Thus there is sood reason to believe that 


our universe does in fact contain a closed trapped Surface. 


it should be pointed. out that the possession of a closed 


trapped surface is a large scale property that does not depend 


on the exact local metric. Thus a model thet had local irrezul- 


arities, rotation and shear but was similar on a large scale 
at the present time to a Robertson-Walker model. would nave a 
closed trapped surface. 

Following Penrose it will be shown that space-time has 
a Singularity if there is a closed trapped surface and 


CE) E y O for any observer with velocity 


(8) there is a global time orientation 


(h) the universal covering space has a non-compact Cauchy 


rr? 
Li 


surface 


PROOF 


Pss 
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Assume space-time is singularity free. Let F be the 

set of points to the past of H^ that can be joined by a 
smooth future directed time-like line to no or its interior 
P2. Let pÍ be the boundary of F'. Local consüeravions show 

ls qus 3 a | T T vm medi Ge 
that D^ - T" is null where it is non-singular and is 
generated by the outgoing family of past directed null geodesic 
segments which have future end-point on no and past end-point 
where or before a singular point of the null geodesic 
congruence. wince at ne , toe convergence, fre? and 

^ 

Since Re iy 2 0 by (f), the convergence must 
become infinite within finite affine distance. "hus 
will be compact being generated by a compact family of compact 


segments. Hence B? will be compact.  Penrose's method is 


; ) ' i 2 4 
then as follows: approximate B^ arbitrarily closely by a 
smooth space-like surface and project B^ onto H^ by the 


normals to this surface. this gives a many-one continuous 


A x 
mapping of 8^ into H^. Since B? is compact, its image B^ 


Z 
must be compact. Let d(Q) be the number of points of B7 


mapped to a point Q of i^. UO) will change only at the 


Z 
intersection of caustics of the normals with H^. Moreover, 


by continuity cL Q) can only change by an even number. 


| ig^) 
But of (73) = 1 since this is the identity and d(H- = Oy, 
This is a contradiction, thus the assumption that space-time 
is non-singuler must be false. An alternative proceedure 
which avoids the slightly questionable step of approximating 
D" by a space-like surface is possible if we adopt condition 
(e) on the nature of the matter, then B? may be projected 
continuously one-to-one onto Hu? by the flow-lines.  -his 


again Leads to a contradiction since pÍ is compact and H 


isa not. 
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-mosphism-of—(]3- into H2, Purther,—by—Ce- 
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6. singularities in 'Closed' Universes 


"m Cram p «m Labo tu. od UT. xd XA و‎ eee 
"nere is a singularity in every model which satisfies 


(a), (g) and (i). 
(i) There exists a compact Cauchy surface g’ whose unit 
Ox 
normal V has positive expansion everywhere on H^. 
vor the proof it is necessary to establish a couple of 
lemmas. Assume that space-time is singularity-free. The 
| following result is quoted without proof, it may reecily be 
derived from lemmas proved in reference 11. 


if o and 9. are conjugate points along a geodesic y 


and. X is a point on 1 not in pa a SUAM have a conjugate 


point in Pa. | 


An immediate corrollary is that if q is the first point 


along y conjugate top and y is in pq then y has no 

conjugate points in pq . Also since tne result that x 

has a conjugate point in pq can only depend on the values 

of Q in pq, any irrotational geođesic congruence including 
Mn 

the geodesic ¥ must have a singular point on y in pa. 

Thus if q is a point on H^ and 1 is the geodesic normal to 

M? through q, then a point conjugate to q along Fi cannot occur 

until after a point conjugate to M^. 

LE M^ is a complete connected space-like surface which 

intersects every time-like and null line from a point p, we 

may define & function over M? as the square of the geodesic 

distance from p which is taken as positive if the geodesic is 

time-like and negative if the geodesic is space-like. e call 

this the world function O with respect to p. For the closed 

set of values © ZO >60 will be & continuous (An 

ceneral multi-valued) function over M^. A time-like Zzeoicsic 
1 from p will be said to be critical if it corresponds 


PE M | 
to à value of G for which y 1 NT d. E 2 § 4, 3) 


PSU 


M 8 
where (© are three independent vectors in M?, Clears 
(7 


a critical geođesic must be orthogonal to M^. A geodesic 


H 
= 


which is critical will be said to be maximal if it corresponds to a 


local maximum of e 
Lemma 1. 
À geodesic Ó cannot be maximal for a smooth M? if there is a 
point X conjugate to M? but no point conjugate to q on f in qp, 
5 


where q is the intersection of A and M^, 


Let f and g be the Jacobi fields along y which vanish at X 


m m 
and f respectively. They may be written 
n 
f = A(s) f/q P 
m mn 
n 
g = B(s)g/a. 
m mn 1 
^ d I i l h j [ o | 5 A 
Then عاك اليس بأ م‎ T- alg / ١ must be positive for 
4 NAD Y (MS L 5 
any h since if it were negative for any h by taking a =-bhh 
mb mn 
A 


beyond q, it would be possible to have a point y on C beyond q 
conjugate to A before a point conjugate to P  . If it were zero 
X would be conjugate to P . This shows that the surface at q 

of constant geodesic distance from pP lies nearer to f in every 
direction than the surface 24 of constant geodesic distance from 
/ 


y does. Since X is conjugate to M^ the surface at q of 
| 2 


P 
constant geodesic distance from [^ lies closer to 2 in some 


2 


x , 
does. Hence Q is not maximal. 


direction than M 


Buta‏ ——,35—3€— هنا 
( 3 
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s 
i If M^ is compact or if the intersection of all time-like 


and null lines with M? is compact, © must have a maximu:: 


value, thus there must be a geodesic normal to M’ through p 
longer than ¥  . We use this to prove another lemas. 
Lemna 2 

If p lies to tne future (past) on a time-like geodesic 


m 


3 through q, beyond a point z conjugate to q, and there 


^ 
exists a compact Cauchy surface i^ through q, then there 


< 


1 must be another time-like geodesic from p to q longer than y 


Let y be the last point conjugate to q on y before p. 
Let x be the nearest point to p conjugate to p in pq. Lect 


r be à point in yx. Let K? be the set of points which have 


a future (past) directed geodesic of length rq from q. then 


NON 


T | | » E 
will be a space-like hypersurface through Fé Let F 


K 
be the set of points which nave at least one future (past) 
directed geodesic from q of length greater than ra. ‘then tne 


Ji 


. F: T" : ; dj. ; 
boundary of F , ye K^. Since p is in F and since every 


past (future) directed time-like and null line from p intersects 


wA NON NEN | " 
which is not in F , they must also intersect J^. Let L^ 


JÍ 2 10 


and these lines. Since dH’ is 


HÎ 


be the intersection of 
compact, L^ must be compact. Consider the function O° with 


5 


respect to p over X^". Its maximum iust lie in the compact 
region x. put, by the previous lemma y is not maximal, 
moreover, local considerations show that a singular point in 
the surface "E cannot be a maximum of O . ‘thus tne 
maximum value of g^ must occur for a geodesic from p ortho- 
sonal to yh This must also be à goedesic from o to q oí 
length greater than Y a 

Usins these two lemmas the theorem may be proved.  oince 


2 


are converging 


2 


everywhere on H^, there must be a point conjugate to H^ a 


the future (past) directed normals to H 


finite distance along each future (past) directed geodesic 
normal. Let $ be the maximum of these distances. Let o 
be à point on a future (past) directed geodesic normal av 


a distance greater than S . Consider the function O with 


respect to p over the compact surface H~. Let be the 


9 


at the point q, wnereO nas its 


maximum. There must be a point conjugate to H^ along A in ap. 


geodesic from p normal to H 


But if there is no point conjugate to q along A in ap, then 
۸ cannot be maximal by the first lemna. If however there 
is a point conjugate to q along À in qp, then tnere must be 
a longer geodesic from q to p by the second lemma. thus 
is not the geodesic of maximum length from H to p. UuUhls.is 
a contradiction which shows that the original assumption tnat 
the space was non-singular must be false. 

his proof could also be used to Snow the occurrence 
of a singularity in a model with a non-compact Cauchy surfece 
provided tnat the expansion of its normals was bounded away 
from zero and provided that the intersection of the Vauchy 


surface with all the time-like and null lines from & point 


was compact. 
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